

%% parameters
load SMM
load('SMM', 'jj')
fprintf('verif poids=1');
p0 = p0_SMM;
mu_c = mu_c_SMM;
%sig_c = sig_c_SMM;
sig_eps_a =jj*1/250;
p_a = p_a_SMM ;
weight_j=0.57;


      [f2, avfracup2,    med2,    interq2, kur2]=geNCalvoPlus_SMM(p0, mu_c, sig_eps_a, p_a, weight_j);
    
    
    
    